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Abstract The stationary phase point (SPP) method in one-dimensional case 
is introduced to treat the diffractive scintillation. From weak scattering, where 
the SPP number N=l, to strong scattering (N3>1), via transitional scattering 
regime (N~2,3), we find that the modulation index of intensity experiences the 
monotonically increasing from to 1 with the scattering strength, characterized 
by the ratio of Fresnel scale rp to diffractive scale r^is- 
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1 INTRODUCTION 

As the radio waves propagate in the interstellar medium (ISM), the diffraction and refraction 
introduced by small-scale and large-scale inhomogeneities lead to the flux variations or scin- 
tillations (review see Rickett, 1977, 1990, 2006; Narayan 1992). The radio wave interference 
fringes are seen evolving with time, due to the motions of source, observer and ISM, which 
can be explained in terms of wave scattering through a random medium, with the electron 
density inhomogeneities being described by a power-law spectrum, close to the Kolmogorov 
spectrum (Armstrong, Rickett and Spangler 1995; Cordes, Weisberg and Boriakoff 1985). As 
for the characteristic scales responsible for the scintillations, Fresnel length scale is defined by 
rp = y/z/k, where z is the distance between the scattering screen and the observer's plane and 
k = 2ir/\ is the wave number, and the diffractive scale rdifF is defined by writing the phase 
structure function in the forms D(r) — {r/r^g) 01 for a Kolmogorov spectrum of turbulence (see 
e.g. Narayan 1992 for a review). In the weak scattering regime, r-p/rss << 1? the flux variation 
is interpreted in terms of weak focussing due to phase curvature on the scattering screen on 
scale ~ rp. Whereas, in the strong scattering regime {r-p/r&g » 1) there are two variation 
scales, small diffractive scale r^g and large refractive scale r ro f = r F /r<jig, caused by interfer- 
ence between the many coherent patches of size r<jiff over the scattering screen. However, in 
the transitional scattering regime, r-p/r^m ~ 1, which is applicable to most intra-day variable 
extragalactic radio sources observed at frequencies between 1 and 10 GHz (see, e.g. Walker & 
Wardle 1998), the flux variation scale is a little fuzzy because rp, rdie and r re f are very similar, 
perhaps a mixing variation scale may be produced. The two time scales of spectrum intensity 
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Fig. 1 The illustration of one dimensional scattering screen. 



from the observations of pulsars reveal variations of minutes to hours as a diffractive scintil- 
lation (see e.g. Scheuer 1968; Manchester and Taylor 1977; Lyne & Smith 2005) and days to 
months as a refractive scintillation (see, e.g. Rickett 1984 for reviews; Kaspi & Stinebring 1992; 
Wang et al. 2005). 

In the transitional scattering regime, where both geometric optics and wave optics are 
equally important, with r^g comparable to rp, the structures on the scattering screen have a 
focal length comparable to the distance from the scattering screen, so the role of the caustic 
focussing will effect here. 

The paper briefly introduces the unified descriptions of modulation index from weak, tran- 
sitional and strong scattering by the stationary phase point method (SPPM), which has been 
previously paid attention by Gwinn et al. (1998) and recently applied to interpret the parabola 
phenomena occurring in the pulsar secondary spectrum (Walker et al. 2004). 



2 THE STATIONARY PHASE APPROXIMATION 

We consider the wavefield from a point source at infinity to be incident up on a single thin, 
one-dimensional scattering screen, as described in figure [TJ For an incident wavefield of unit 
amplitude the wavefield observed at the central position on the line a distance z (rp = z/k) 
from the scattering screen is (e.g. Born & Wolf 1980) 



u(z) 



1 



dxexp 



2r| 



To solve this we employ the method of stationary phase (e.g. Mandel & Wolf, P. 128), which 
states that an integral of the form, 



F(k) 



f{x)e lk9{x) dx, 



(2) 



has an approximation containing contributions from critical points inside the boundary of in- 
tegration (—00,00): 



F(k) 



1/2 



(3) 
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where ej = e ±I7r / 4 is according as g"{xj) > (< 0) for + (-) and Xj are the points of 
stationary phase which satisfy the condition g'(xj) = 0. Applying the stationary phase 
approximation and g(x) = k^ 1 ((j)(x) +x 2 /2r|) to Eq.fll]), one has 

N 

u{z) tt^Uj, Uj = \ = exp[irl(f>' 2 (x J )/2 + i(j>{xj)]. (4) 

Jit v|1+»f^(*j)I 

It is noted that Eq.(3) or Eq.(4) is an approximate result, and it is assumed 
that g" (x) does not go to zero. The independent variables are the phasor cf>, its first 
derivative <fi', and the total number of them is 2N, where N is the SPP number. 

The mean intensity of the unit amplitude wave, < / >= 1, is identified as the average value of 
the second moment of the amplitude < uu* >. Adding a subscript N to denote the number of 
SPPs, in the SPP approximation this is 

N N N 

< i > N = ^ < Uiu *j >= < UiU * i > + ^2 < Uiu *] >= l - ( 5 ) 

i,j=l i=l i^j=l 

Thus the number of SPP is the maximum integer of N solved from Eq. ^ . While considering 
a statistical model in which each SPP is assumed to come from a distribution, all SPPs are 
statistically equivalent. Then, for example, < UiU* > may be written as < u\u\ >, and the sum 
over i — l,...N gives N < U\u* >. Only two independent averages appear in the intensity, 
and we write these as 

U u =< uiul >, U 12 = (< uxu* 2 > + < u 2 ul >)/2, (6) 
so that Eq.© reduces to 

< I > N = NU n + N(N -1)U 12 = 1. (7) 

Similar to the procedure of treating the intensity, the intensity square is obtained as follows, 

< I 2 >jv= NUuu + 4N(N - l)t7 m2 + N(N - l)(U 1122 + U 1221 + U 1212 ) 

+N(N - 1)(N - 2)(4C/ 1123 + 2U 1213 ) + (N - 1)(N - 2)(N - 3)U 123i , (8) 

where the terms, such as Uuu and Uu 22 , are the averaged values of the products of amplitudes 
and generally denoted by, 

Ujkim =< jklm >=< u.jU* k uiu* m >, (9) 

which satisfies, 

< jklm >*—< kjml >—< mlkj >, (10) 

so we have Uu 22 = U\ 22 \ = UuU 22 = U\ x . Henceforth, the modulation index m of diffraction 
is defined by 

m = v / < I 2 >n -1 • (11) 
2.1 The calculations of the intensity in SPP 

The power spectrum of the phase inhomogeneities, <&(&), is assumed to have the following form 
(see e.g. Blandford & Narayan 1985; Goodman & Narayan 1985) 

L , k < 1/L 

$(fc) = { L (Lk)- p , l/L < k < l/l , (12) 
0, k>l/l 
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where (3 = 8/3 is an index for a Kolmogorov spectrum of turbulence with the structure constant 
L and I (L) is the inner (outer) scale. In addition, we assume the phase parameters to follow 
the Gaussian distribution P(ip)i where i/j represents the phase variables <fi, <fi' and <f>" , 



i 



2no* 



: exp 



(13) 



where is a standard deviation of the parameter ip, calculated from the spectrum function of 
Eq. (fT2"]) (see Melrose & Watson 2006). Furthermore, the terms in intensity and intensity square 
expressions of Eq. ([7]) and Eq. ([5]) can be calculated by an integral of the probability distribution 
over the 6 random variables (see Melrose & Watson 2005), for instance, 



< jk >= (u jU * k ) = I [ 





171-/21 



do 

P( 



1/2 



l + r^| 



1/2 



) 

exp [iZjk] 



do 



ill iTT/2 



with 



■'jk 



P^^k, (ft) = P(^)P(^)P(<^)P(^)P(^0«) 



(14) 

(15) 
(16) 



3 RESULTS AND DISCUSSIONS 

The numerical results of the modulation index versus the scattering strength rp /ram are shown 
in figure [2j We find that the modulation index increases with the scattering strength continually 
from the weak scattering regime, via the transitional scattering regime, to the strong scattering 
regime where the modulation index is fully saturated. The influence by the inner scale on 
the modulation index is studied, and we find that the increasing of inner scale makes the 
modulation index decrease a little at the transitional scattering regime but little effect in the 
strong scattering regime. The trends of the modulation index in the weak and strong scattering 
regimes are similar to those obtained by the approximated treatments (see. e.g. Narayan 1992; 
Rickett 1977, 1990), however the SPP method provides the descriptions of modulation index at 
the transitional scattering. If N=l, a weak scattering case, we have 

m « [/miM 2 ! - 1 , (17) 

from which we can obtain that the modulation index increases with the scattering strength 
Tf/t^s, similar to the case discussed by Salpeter (1967) (Melrose & Watson 2005). In the 
transitional scattering the number N is about 2 or 3 but depends on the choice of inner scale. 
If N — > oo, a strong scattering case, we find that the interchange terms of amplitude products, 
for instance < jklm > with j,k,l and m being not all same, are very small and Umi ~ t/1122 ~ 
E/1221) ~ 1/N 2 . Therefore, the following approximation is preserved, 

m =< I 2 > N -1 ~ NUnn + N(N - 1){U 1122 + U 122 i) -1^1, (18) 

a fully modulated scintillation, which has been obtained and discussed by Gwinn et al. (1998). 

Moreover, there exists a discrepancy between our results and those earlier re- 
sults by Goodman and Narayan (2005) who found the modulation index m>l at 
transitional scattering regime, which calls into question the reliability of the SPP 
method. Or it may possibly be due to the 1-D nature of our computation, which 
could be resolved by applying the method to a 2-D screen in a future work. 
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Fig. 2 Modulation index v.s. scattering strength rp/rais- The influence by the inner 
scale of ISM is plotted as indicated in the figure. 
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